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It is shown that the 99% confidence limits from the analyses of the 
data of cosmological and neutrino experiments imply a small marginally 
^ ' allowed region in the space of the neutrino oscillation parameters of 3+1 

. four-neutrino mixing schemes. This region can be confirmed or falsified 

f-^ , by experiments in the near future. 

CN ! The impressive results of the first year of observations of the Wilkinson Microwave 

^ I Anisotropy Probe (WMAP) [1] combined with other recent high precision cosmological 
O I data sets (2dF Galaxy Redshift Survey, Lyman a forest data, CBI and ACBAR microwave 
^ ' background data; see references in Ref. [2]) allowed to derive tight constraints on the 
values of cosmological parameters [2]. In particular, the combined fit of cosmological 
Qh! data performed in Ref. [2] imply the tight bound 

Quh"^ < 7.6 X 10"^ (95% confidence limit) (1) 

^ ■ for the contribution of massive neutrinos to the energy density of the universe {Qu is the 
neutrino density relative to the critical density, h is the Hubble constant in units of 100 
km/s/Mpc), whose dependence on the values of the neutrino masses rrii is given by [3,4] 



93.5 eV ^ ' 

Strong evidences in favor of neutrino oscillations have been obtained in atmospheric 
[5-7], solar [8-13], long-baseline reactor [14] and accelerator [15] neutrino experiments. 
The results of all these neutrino experiments can be explained in the framework of stan- 
dard three-neutrino mixing (see Refs. [16,17]) with the two squared-mass differences 

Am,^„i ~ 7 X 10"^ eV^ , Am^^^ ~ 2.5 x lO'^ eV^ , (3) 

which generate, respectively, solar and atmospheric neutrino oscillations. In this case, 
the bound (1) can be saturated only if the three neutrino masses are almost degenerate 
(see Ref. [18]), and each neutrino mass is bounded by [2] 

mi<0.23eV (95% confidence limit) . (4) 
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Figure 1: Schematic view of the masses in the two 3 + 1 four- neutrino schemes in which 
the isolated neutrino is the heavier one. 

Neutrino oscillations in the framework of three-neutrino mixing cannot explain simul- 
taneously the results of solar, atmospheric, long-baseline neutrino experiments and the 
evidence in favor of short-baseline — > z/g transitions found in the LSND experiment [19], 
because an additional large squared-mass difference. 
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is needed. Taking into account only the need of three independent squared-mass differ- 
ences, the minimal framework that may be able to explain all data with neutrino oscil- 
lations is four-neutrino mixing (see Refs. [16,17]). There are two types of schemes that 
provide the needed hierarchy of squared-mass differences in Eqs. (3) and (5): 2+2 schemes 
with two pairs of almost degenerate massive neutrinos separated by the so-called "LSND 
gap" of the order of 1 eV, and 3 + 1 schemes with a triplet of almost degenerate massive 
neutrinos and an isolated massive neutrino separated by the LSND gap. In spite of the 
apparent large freedom in four-neutrino mixing schemes, which have three independent 
squared-mass differences and six mixing angles, the present data of neutrino experiments 
are so rich that four-neutrino schemes are tightly constrained. In Refs. [20, 21] it has 
been shown that 2 + 2 schemes are strongly disfavored by data and 3 + 1 schemes fail to 
provide an acceptable fit of the data, except for a few marginally acceptable regions in 
parameter space. 

In a very interesting paper [22] Pierce and Murayama considered the marginally al- 
lowed 3 + 1 schemes shown schematically in Fig. 1, in which the isolated neutrino is the 
heavier one, with mass 7714. In this case, only the heaviest massive neutrino saturates 
the cosmological bound in Eq. (1), leading to the limit 



777,4 < 0.71 eV (95% confidence limit) 



(6) 



which straightforwardly implies 



A<SND < 0.50 eV' 



(95% confidence limit) . 



(7) 
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Figure 2: Figure taken from Ref. [21], on which we have superimposed the 95% confidence 
hmit in Eq. (7) (dashed horizontal fine) and the 99% confidence hmit in Eq. (11) (sohd 
horizontal line). 

In Fig. 2 we show this limit (dashed horizontal line) superimposed to Fig. 4 of Ref. [21], 
which shows the LSND-allowed regions at 90% and 99% CL. in the plane of the mixing 
parameters sin^26'LSND, ^"^lsnd) compared with the upper bound for sin^ 2^lsnd from 
short-baseline, atmospheric and solar neutrino data. These bounds are due to the fact 
that in 3 + 1 four-neutrino mixing schemes the oscillation amplitude sin^ 2^lsnd in the 
LSND experiment (as well as in other short-baseline t'^ — > and experiments) 
is given by 

sin2 2^LSND = 4|t/e4nf/M4l'> (8) 

where U is the unitary mixing matrix (see Ref. [16]). The values of |f/e4p and |?7;i4p 
are constrained by the negative results of short-baseline reactor v^. — > and accelerator 
^ii disappearance experiments, whose oscillation amplitudes are given, respectively, 

by 

sin2 2ep,^p, = 4|f/e4|' (1 - |[/e4r) , siu^ 20,^^,^ = A\\J ^^\^ (1 - 1^/^41') . (9) 

Large values of |f//j4p are forbidden by the results of atmospheric neutrino experiments, 
because in order to have oscillations due to Am^tm the muon neutrino must have a 
large mixing with the three light massive neutrinos, which by unitarity implies a small 
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mixing with the heavy isolated massive neutrino. Similarly, large values of |C/e4p are 
forbidden by the results of solar neutrino experiments. The resulting upper bounds on 
sin^ 2^LSND calculated in Ref. [21] at 95% and 99% C.L. are represented by the wiggling 
lines in Fig. 2. 

One can see from Fig. 2 that the cosmological upper bound (7) excludes all the regions 
which are marginally allowed by solar, atmospheric and terrestrial neutrino experiments. 
From this comparison Pierce and Murayama [22] concluded that a four-neutrino mixing 
explanation of the LSND results is strongly disfavored and "the only way to reconcile 
LSND with the cosmological data is to have CPT violation" [23, 24] . 

Here we would like to point out that, although there is a clear tension between dif- 
ferent data interpreted in the framework of four-neutrino mixing, before considering this 
possibility ruled out we should apply more caution. A way to avoid premature conclu- 
sions is to consider experimental limits with high confidence level or high probability, 
which is especially recommended for the comparison of results of different experiments. 
In this spirit, the 95% confidence of the limit in Eq. (7) does not appear to be sufficient 
to rule out four-neutrino mixing. 

Unfortunately, the authors Ref. [2] did not give the values of more robust confidence 
limits. However, it is possible to extract additional information from Fig. 14 of their 
paper, which shows the cumulative probability of Q^h^. One can see that the cumulative 
probability of ^^h'^ is rather fiat above the limit in Eq. (1) and requiring 99% confidence 
may raise significantly the bound. 

The precise values of the cumulative probability of Q^h'^ extracted from the Postscript 
file of Fig. 14 of Ref. [2] are available at ftp://wftp.to.infn.it/pub/giunti/wmap-0302209/ 
wmap-0302209-fl4.txt. The plot in Fig. 14 of Ref. [2] ends at Q^h"^ = 0.010, where we 
find a cumulative probability of 1.4 x 10~^ which correspond to a 98.6% confidence limit. 
For simplicity, let us approximate it with a 99% confidence limit. Hence, in the 3-1-1 
schemes in Fig. 1 we have the bound 

m4<0.94eV (99% confidence limit) , (10) 

and 

^"^LSND < 0.87 eV^ (99% confidence limit) . (11) 

This bound is represented by the solid horizontal line in Fig. 2. 

From Fig. 2 one can see that the three 99% confidence regions obtained from the LSND 
results, from the results of solar, atmospheric and terrestrial neutrino experiments, and 
from cosmological data overlap on a tiny region in the sin^ 2^lsnd~'^'^lsnd plane at^ 

sin^ 2^LSND ^ 2.0 x 10"^ , Am^gND - 0-86 eV^ . (12) 

This region corresponds to the region R2 found in Ref. [26]. 

Although we have to admit that this region is only marginally allowed and its prob- 
ability is rather low, we think that it is still not excluded and may constitute the "last 
CPT-invariant hope for LSND neutrino oscillations"^. 

^See also the interesting Addendum 4 in Version 4 of the hep-ph version of Ref. [25], which appeared 
after completion of this paper. 

^We do not consider here the constraints imposed by Big-Bang Nucleosynthesis, whose status is still 
controversial (see Refs. [27-30]). 
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The predictions of the values of neutrino oscillation parameters in Eq. (12) for future 

experiments arc rather strong. Of course, such region should be measured in the Mini- 
BooNE experiment [31], which is aimed to check the LSND result. Since this region is 
just below the upper bound for sin^ 26'lsnd from short-baseline, atmospheric and solar 
neutrino data, for a value of AmLSNo where such bound is due to the combined results 
of the Bugey reactor Pg disappearance experiment [32] and the CDHS accelerator dis- 
appearance experiment [33], future experiments should find v^. and disappearance at 
the borders of the Bugey and CDHS bounds at Am^ ~ 0.86 eV^: 

sin^ 2^p^^p^ ~ 7 X 10-2 , sin^ 2d,^^,^ ~ 0.1 , (13) 

which correspond to 

|C/e4|2 ~ 2 X 10-^ , |C/^4|^ ^ 3 X 10-2 (14) 

(see also Figs. 3-6 of Ref. [26]). Furthermore, if there are no fine-tuned cancellations 
among the contributions to the effective Major ana mass in neutrinoless double-/? decay 
[34] , this effective mass could be as high as 

|C/e4|^m4 ~ 2 X lO-^eV, (15) 

and could be observed in future experiments [35-38]. The contribution of the heavy 
neutrino z/4 to the effective neutrino mass in tritium beta decay (see Ref. [18]) is 

|?7e4|m4 ~ O.leV, (16) 

which is not much below the 0.35 eV expected sensitivity of the KATRIN experiment [39] 
and may be observed in the future. 

In conclusion, we would like to remark that, although the alternative possibility of 
large CPT violations in the neutrino sector is undoubtedly very interesting, great caution 
should be exercised before ruling out four-neutrino mixing schemes that may allow an 
appealing explanation of all neutrino experiments in the framework of local quantum field 
theory with the standard three flavor neutrinos and an additional very intriguing sterile 
neutrino (see Ref. [40]). Considering 99% confidence limits from the analyses of the data 
of cosmological and neutrino experiments we found a small marginally allowed region for 
the oscillation parameters in 3 -|- 1 four-neutrino mixing schemes, which can be confirmed 
or falsified by experiments in the near future. 
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